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Abstract 



We have investigated edge modes of different multipolarity sustained by 
quantum antidots at zero magnetic field. The ground state of the antidot is 
described within a local density functional formalism. Two sum rules, which 
are exact within this formalism, have been derived and used to evaluate the 
energy of edge collective modes as a function of the surface density and the 
size of the antidot. 
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Typeset using REVTeX 



With the progress of microstructure technology, the study of the two dimensional electron 
gas (2DEG) has evolved to that of laterally confined superlattices made of either electron 
islands (dots) or holes surrounded by electrons (antidots). Much effort has been devoted in 
the past to the study of quantum dots, as compared to that put in the study of quantum 
antidots. One of the goals of their study has been to determine the far-infrared response of 
these semiconductor microstructures, and the formation of compressible and incompressible 
states when a magnetic field B is perpendicurlaly applied. In the case of antidots, which 
is the subject of the present work, experimental evidence of collective excitations sustained 
by these structures has been presented in Refs. A theoretical description based on 

magnetoplasmons in two-dimensional antidot structures has been given in Ref. |4[] which 
compares well with the experimental data of Ref. |2j . Recently, the existence of compress- 
ible and incompressible strips at the edge of antidots has been determined by far-infrared 
spectroscopy 0. 

We have started a systematic study of the structure and collective far-infrared response 
of antidots, whose aim is to achieve a level of sophistication in the description of these 
systems similar to that attained for quantum dots. As a first step, we present here results 
at zero magnetic field obtained within Density Functional Theory (DFT). To some extend, 
the present study is similar in scope to that carried out in Ref. on the surface excitations 
of cavities in 3D metals. The B ^ case, which requires a rather different and far more 
complex approach, will be presented in a forthcoming paper. 

We have modelled an antidot of radius R in a 2DEG of surface density n s by a positive 
jellium background of density n s Q(r — R) to which we have added for the reason that will be 
given below, a parabolic potential barrier of the type V p (r) = m*uj 2 {R 2 — r 2 )/2 acting on the 
electrons for r < R. We shall call V ex t the sum of the jellium and V p potentials. The ground 
state (gs) of the antidot is obtained solving the Kohn-Sham (KS) equations as indicated 
for example in Ref. |7|]. The problem is simplified by the imposed circular symmetry, and 
only the radial KS equations have to be considered to determine the electronic radial wave 
functions RAr): 
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where m* is the electron effective mass which together with a dielectric constant e are 
characteristics of the semiconductor. For example, e = 12.4 and m* = 0.067 m e in GaAs, 
which we have chosen for the numerical applications in view of the existing experimental 
data of Refs. [@||. I = 0, ±1, ±2... is the angular momentum about the z axis perpendicular 
to the plane of the antidot, and r is the radial variable in that plane. The single electron 
potential V(r) is made of V ex t, of the Hartree electron-electron potential and of the exchange- 
correlation potential. The correlation potential has been obtained from the correlation 
energy of Ref. || in a local density approximation. Some times we shall use effective atomic 
units. In this system of units, the length unit is the Bohr radius times e/m*, and the energy 
unit is the Hartree times m*/e 2 , denoted here as Oq and H*, respectively, e is the dielectric 
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constant, and m = m*m e is the electron effective mass. For GaAs we have e=12.4 and 
m*=0.067. Consequently, a* ~ 97.9A and H* ~ 11.9 meV ~ 95.6 cm -1 . 

Physically acceptable solutions to Eq. ([[J) have to be regular at r = and behave 
asymptotically as Ri(r,k) ~ Ji(kr) + tan(5i)Yi(kr) , where Ji and Y\ are the integer Bessel 



functions of first and second kind |J, and k = y2m*(E — V^/Ti 2 . The wave number k has 
to be k < kp = y/27rn s . Taking into consideration the spin degeneracy, the electron density 
p(r) is obtained as 

M=T^J^ r Rl(r,k) d k. (2) 

We have checked that the number of points used to carry out a Gaussian integration over k, 
and the maximum \l\ employed in Eq. (Q) lead to stable results. Typically, some 1000-1500 
wave functions have been computed to obtain a density. 

We show in Fig. (|]) the electronic densities corresponding to antidots of R = 7.5 Oq and n s 
= 0.05, 0.1, 0.2, 0.3 and 0.4 (dp) -2 . They span the size and density range of those fabricated 
and experimentally studied in Ref. [Q, and have been obtained using u = 0.09H*. Figure 
(|) shows the V(r) KS potential corresponding to n s = 0.1 and 0.4 (ag) -2 - Also indicated 
are the Fermi energy and the contribution of V p (r) to V(r). 

Looking at the position of the Fermi energy with respect of the top of V(r), it is not 
surprising that a potential barrier of a kind or another is needed to prevent the electrons from 
spilling in the antidot too much, producing an unphysical representation of the experimental 
device. Actually, we have found that if u is set to zero, high density antidots would have 
non-zero electron densities even at r = 0. This is illustrated in Fig. (|3]) for the R = 180 
nm, n s = 9 x 10 11 cm -2 antidot of Ref. p. In this figure, the dashed line density has been 
obtained setting uj q = 0. 

These results can be employed to determine the B = far-infrared multipole response of 
antidots. To do so, we rely on the formalism described in detail in Ref. []TU| . For the present 



purpose, it consists in obtaining the so-called mi and 7713 sum rules (SR) for an excitation 



operator Ql- We have that |fTT|: 



MQl) = \(o\[Q£,[h,Ql]]\o) 



m 3 (Q L ) = ~(0\[[H, [H,Q+]],[H,Q L 



(3) 



where |0) is the gs of the system. These SR have been extensively studied in the literature 



rT,12|. As indicated in these references, if only a collective state is contributing to the 
strenght function, Es(Ql) = (%/wii) 1 '' 2 represents the average excitation energy. This is 
the situation experimentally found for antidots at zero magnetic field. 
The operator Ql is taken to be 



N j 



Ql = T,-l^ ■ (4) 
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This choice is inspired in that (qr)~ L e lLd is the small q expansion of the function YL(qr)e tLe , 
which is the restriction to the z = plane of the irregular solution of the Laplace equation in 
cylindrical coordinates. Following Ref. ||10|| , a lengthy but straightforward calculation yields: 
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where 



m 3 (Q L ) = m 3 (T) + m 3 (ee) + m 3 (V ext e) 



m 3 (T) = 2nL 2 (L + 1) j dr [Lr(r) + 2(L + 2)A(r)] 
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The definition of the densities r(r) and A(r) and of the function El can be found in Ref. |10 



and the primes on the density denote r-derivatives. The jellium potential V + (r) entering in 
Vext is: 



Vjr) 
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r < R 



R^ E(r/i? 00 ) - r E(R/r) + r (1 - (i?/r) 2 ) K(R/r) r> R. 
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In the above equations, K and E are the complete elliptic integrals of first and second kind, 
respectively flI3| , and R^ represents a large r value. In practice, it is the largest r used 
in the structure calculation, which we have also taken as the point where the asymptotic 
behavior of Ri(r) holds. We want to point out that the two Coulomb diverging terms in 
m 3 (ee) and m 3 (V ext e) cancel each other. 

The present formalism can be applied to antidots with the restriction that p(r) vanishes 
in a small disk around the origin. In practice, this is not a limitation, as can be inferred 
from Figs. (see also Fig. H). Some technical details about how the above integrals are 

handle can be found in Sec. IV of Ref. . 

For a large antidot, the electronic density is constant everywhere apart from a narrow 
region along the border of the hole. Following the method outlined in Ref. 0] , it is easy to 
show that E$ yields the classical hydrodynamics dispersion relation for edge waves, namelly 
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Es = Lo(q) ~ y2n s q ln(g /<?)- It is also worth to notice that the induced (or transition) 
density associated to the operator Ql has the form JlO 



p tr (r) ^ L-^ p\r) (11) 



that manifests the edge character of the excitation. 

Figure (|j) shows the frequency of the L = 1 mode as a function of n s compared with the 
experimental points of 0. For completeness, we have also plotted the results corresponding 
to L = 2. We have checked that similar results are obtained using as potential barrier the 
parabola V p (r) = m*ujQ(R — r) 2 /2 for r < R with ujq = 1H*. One can see that the agreement 
with experiment is good. Furthermore, our calculation yields a frequency of ~ 68 cm^ 1 for 
the R = 180 nm, n s = 9 x 10 11 cm~ 2 antidot, in good agreement with the experimental 
findings of Ref . . 

We have also studied the size dependence of the mode energy. Fig. (||) represents the 
electronic densities for antidots of R = 10, 15 and 20 a^, and a surface density n s = 0.2 
(czo) -2 . The frequencies of the L = 1 and 2 modes are shown in Fig. (|6]) as a function of 
1 / \fR. They exhibit a distinct R dependence, indicating a clear departure from parabolicity 
of the confining potential, i.e., a physical situation where the generalized Kohn theorem does 
not apply. 

If the electronic density p(r) is approximated by a quasi-step function, an analytical 
expression can be obtained for E\. Proceeding as in Ref. fl(| one gets: 



.2 L(L + 1) A L 



l ln (j?) + l-£ 



(12) 
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where a represents the width of the quasi-step function, and the precise value of 7 depends 
on the way the electronic density goes to zero [14]. This equation tells one that the frequen- 



cies have a 1/yR linear dependence if the Coulomb energy contribution dominates. For 
reasonable values of 7/a, it happens for any realistic value of R. 



Equation (|TJ) is too crude an approximation because of the quasi-step function density 
used to get it. Much better results can be obtained if Eqs. (|5|-|9]) are evaluated in the 
Thomas-Fermi approximation, i.e. r(r) = 7ip 2 (r), A(r) = r(r)/2 flOH , after having fitted the 
KS density to a generalized Fermi function of the kind 



1 _|_ e (r-R)/a 

For n s = 0.2 (ag) -2 ; the KS densities are well reproduced on average taking v = 1.1 and a 
= 0.38 <2q (see Fig. ||). The mode frequencies are represented by the lines in Fig. (|]). 

In conclusion, we have shown that Density Functional Theory is able to reproduce the 
zero magnetic field frequency of antidot edge modes in quite a similar way as it does for 
quantum dots. Although a satisfactory description of the collective spectrum of antidots 
can be achieved using a magnetoplasmon approach H, an alternative method based on a 
more microscopic approach such as DFT is needed to discuss other interesting problems, 
such as edge reconstruction and the formation of compressible and incompressible strips at 
the antidot edge ||. 
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FIGURES 



FIG. 1. Electronic densities as a function of r for antidots of R = 7.5 ag and n s = 0.05, 0.1, 
0.2, 0.3 and 0.4 (ag) -2 . Also shown are the jellium densities (dotted lines). 

FIG. 2. KS potential V(r) as a function of r for antidots of R = 7.5 ag and n s = 0.1 and 0.4 
(ag) -2 (solid lines). The dashed lines represent V(r) without the V p (r) contribution. The horizontal 
solid lines indicate the Fermi level, and the vertical dotted line, the radius of the antidot. 

FIG. 3. Same as Fig. 1 for R = 18.35 ag and n s = 0.86 (ag) -2 . The dashed line density has 
been obtained setting ujq = 0. 

FIG. 4. Mode frequency as a function of the electron surface density for L = 1 and 2. The 
crosses are experimental data from Ref. [2], and the lines are drawn to guide the eye. 

FIG. 5. Electronic densities as a function of r for antidots of R = 10, 15 and 20 ag, and n s = 
0.2 (ag) -2 . Also shown are the jellium densities (dotted lines). The thin solid line represents the 
parametrized density, Eq. [13], for R = 20 ag. 

FIG. 6. L = 1 and 2 mode frequencies as a function of R" 1 ! 2 for n s = 0.2 (ag) -2 . From right to 
left, the dots correspond to R = 7.5, 10, 12.5, 15, 17.5, and 20 ag. The lines represent the results 
obtained using the density Eq. [13]. 
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